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V. — On a Classification of Elastic Media, and the Laws of Plane Waves propa- 
gated through them. By the Rev. Samuel Haughton, Fellow and Tutor of 
Trinity College, Dublin. 



Bead January 8, 1849. 



IN a Memoir on the equilibrium and motion of solid and fluid bodies, pre- 
sented to this Academy in May, 1846, 1 deduced the laws of such bodies from 
the hypothesis of attracting and repelling molecules ; since that time I have 
been led to consider the general laws of continuous bodies, without making 
any such restriction as to the nature of the molecular action. The present 
paper contains the results I have arrived at in this investigation, and may, per- 
haps, be considered interesting on account of the classification suggested as ap- 
plicable to all elastic media. It consists of five sections ; the first contains the 
general equations applicable to all media, and the properties of plane waves 
transmitted through them, which are readily deduced from an extension 
which I have given to a theorem originally stated by M. Catjchy, for a parti- 
cular case ; the second, third, and fourth sections contain respectively the laws 
of the three groups into which elastic media may be divided ; these three groups 
consisting of, — first, bodies whose molecular action consists of exclusively nor- 
mal pressures ; secondly, bodies whose molecular action produces exclusively 
tangential forces ; thirdly, bodies composed of attracting and repelling molecules. 
The fifth section contains a comparison of the mechanical theories of light pro- 
posed by Mr. Gbben and Professor Mac Cullagh, with some observations on 
the present state of the science of physical optics. Whatever theoretic objec- 
tions may be made to the application of the theory of elastic media to optics, 
none such exist as to its application to solid and fluid bodies The mathematical 
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investigations which, in the case of light, must be hypothetical, are, in the case 
of solids and fluids, essentially positive, and may be made the object of direct 
experiment. A general inquiry into the laws of elastic media, is an interesting 
application of rational mechanics, and although it must necessarily include cases 
purely hypothetical, it is not, therefore, to be considered as unimportant. In 
this respect it is analogous to an inquiry into the general theory of central forces, 
the importance of which is not confined to the investigation of those laws, of 
which examples occur in nature ; these are undoubtedly the most important, 
but the theory of central forces, considered as a branch of mechanics, would 
be incomplete, unless extended to all possible laws of central force. 



SECTION I. — GENERAL EQUATIONS. 

The formula of virtual velocities, which contains, as shown by Lagbange, 
the conditions necessary to be fulfilled in the interior and at the boundaries of 
a continuous body* is the proper starting point for a deductive theory of the 
mechanical structure of bodies. Every hypothesis which may be made, and 
every fact which experience has discovered, respecting the molecular consti- 
tution of bodies, may be expressed in its most simple form by the aid of this 
formula ; which, by its flexibility, and the facility it affords for deducing theo- 
retic results, becomes of more importance in questions of this nature than in 
other mechanical problems. In order to express by means of it the condi- 
tions of equilibrium of a continuous body, it is necessary to distinguish the 
forces acting at each point into two classes, molecular and external forces ; in- 
cluding among the external forces the resultants of the attractions of the other 
points of the body, since these attractions arise from gravitation, and must not 
be confounded with the molecular forces. The formula of virtual velocities 
must also be stated in such a manner as not to involve any hypothesis as to 
the nature of molecular forces, so as to possess the requisite degree of generality. 
If the problem be dynamical, we must then add accelerating forces equal and 
opposite to those actually employed, so as to destroy the motion at each point, 
and consider the problem as one of equilibrium of forces. These negative accele- 
rating forces must be considered as external forces. 
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The forces being thus divided into two classes, the formula will consist of 

two parts, 

2 (Pip + Pip' + &c. ) + 2 ( Qlg + Qlg' + &c.) = 0. 

P, P, &c, denoting the external forces, and Q, Q', &c, denoting the molecular 
forces. If (x, y, z) denote the co-ordinates of a point at any instant of the mo- 
tion, and a virtual displacement (8ar, 8y, Iz) be conceived, these quantities will 
be functions of four independent variables, which will be the initial values of 
(x, y, z), and the time. Let («', i/, «/) denote the accelerating forces ; hence 

2 (Pip + P V + &c.) = Iff { (X - u 7 ) Ix + ( Y- v') ly + (Z- w') Iz j dm. 

I shall assume that the virtual moments of the molecular forces depend upon 
the differential coefficients of (8a?, ly, Iz), by means of the following linear 
equation,* 

<*+«*+* -P. £+*$+*£ 

Inserting these values into the equation of virtual velocities, and integrating by 
parts, we find 

= JJJ{/» (X -u')lx + p( F- v') ly + p(Z- w') Iz \ dxdydz 

+ \\( P x lx + Q,8y + R,lz) dydz, ( 1 ) 

+ JJ(P 2 8a: + Q 2 ly + B 2 lz) dxdz, 
+ fl(P»8a? + Qfiy + R£z) dxdy. 

The equations of motion being determined by the triple integrals, and the con- 
ditions at the limits by the double integrals. The equations of motion formed 
from this equation will be 

* The ground of this assumption is the fact, that molecular forces depend upon the relative dis- 
placements of the particles, and not on their absolute displacements. 

o2 
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xr-»o=§ + f + §, (») 

^ ' dx dy dz ' 

These equations are the same as those deduced by a totally different method 
(which will be given presently), and involve no restriction as to the extent of 
the deviation of (#, y, z) from their original values (a, b, c). 

If we suppose x = a + £, t/ = b + ij, z = c + £, the expression for the mo- 
ments of the molecular forces will become 

Q$q + QV + &c. = PM + P&h + Pago, 
+ Qffli + Q 2 8ft + QM 
+ Rfyi + R 2 hy 2 + .ffj&y, ; 
where 









a- dr > 
^~Tx' 




* ds' 






k 

*=£• 



If we restrict the molecular forces by the condition 
we shall have the relations 



da t ' dai ' 




a- dV o- dV 


^"tfft' 


Ri = d Z, R t JJ- 

dy^ * dy, ' 


«7s 



F=F(ai, as, a,, ft, ft, ft, 71, y 2 , ?$)• 
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The kind of motion -which it is the object of this paper to investigate is of 
the kind commonly called small oscillations ; and for this kind of motion it is 
not necessary to use the most general equations, or to consider the unknown 
quantities of the problem as functions of (x, y, z, t). We may use, instead of 
(x, y, z), the co-ordinates (a, b, c) of the position of rest of the molecules. In 
fact, any differential coefficient of a function <p, taken with respect to (a, b, c), 
will be expressed by the equation, 

d<f> _ d<f> dx d<p dy d<f> dz 
da~ dx da dy da dz da ' 

but, since x = a + £, y = 5 + rj, z = c + £, we obtain 

dx _ <Z| dy _drj dz _d£ 
da da' da~ da' da da ' 

Hence, neglecting quantities of the second order, we find 

d<f> __ d<j) 
da dx' 

and similarly for the other differential coefficients. 

In the remaining part of this memoir (unless the contrary be expressed), I 
shall, therefore, consider (x, y, z) as the co-ordinates of the position of equili- 
brium of the molecules, and (£, »/, £) as the small displacements of the molecule ; 
the element of the mass will be expressed by the equation dm — edxdydz, 
where e denotes the density, not considered as a function of the time, since 
(dxdydz) denotes the original element of the volume. 

Two kinds of waves can pass through such a body as water ; one, a surface 
wave, depending on the action of gravity for its propagation ; the other, such 
a wave as propagates sound, and does not directly depend on external forces. 
This latter is the kind of wave described in this paper. The equations pecu- 
liar to it will be found by omitting (X, Y, Z) from the general equations ; but 
though the external forces are not explicit in the formulae, yet as they affect 
the density and structure of the body differently at different points, though they 
do not directly affect the wave, we must introduce them implicitly by rendering 
e and the coefficients of V functions of (X, Y, Z), and therefore of (x, y, z) 
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The equation of virtual velocities thus modified will become, considering 
(x, y, z) as the positions of rest of the molecules, 

l\^^^ l ^W^)^dydz = ^Vdxdydz. (3) 

As V is a function of the quantities (oi, aj, a s , &c.) of a given form, we shall 

have 

cT _ dV e dV^ dV c 
8 V = j— fea, + -j— 8a2 + j — fcos 
Oai <zo2 acts 

<* F ™ <*F*„ <**%„ 

dV ,dV dV 

dyi dy 3 2 dy 3 

Substituting this value in equation (3), and integrating by parts, we obtain 
f 6 (S^ + ^ ^ + M K)d*dydz = W>Vdxdydz 

_ fff/A ^Z. £ lE+£ . ^T\ s ^^ & 

JJJ \efo; ^o t dy dot dz da 3 ) 

m(d dV d dV d dV\^ , , . 

rrr/rf dV d dV d dV\^ . . , 
-^{T,-W l + ry'dy- i + Tz-dy-J^ dxd ^' 

The double integrals, as usual, denote the conditions at the bounding surface, 
and the triple integrals give the general equations of motion in the interior. 
The laws of propagation depend upon the triple integrals, and the laws of re- 
flexion and refraction depend upon the double integrals. 



The equations 


of motion derived from (4) will be, 
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dV 

da n ' 

dfh ' (5) 

dV 

dy 3 ' 

These equations contain the laws of propagation of every variety of wave 
not depending on external forces ; and as the function has not been restricted 
by any hypothesis as to the law of molecular action," they will be the dynamical 
equations of propagation of sound in air, water, solids, and of light, if we adopt 
the undulatory hypothesis. In all these cases, the difficulty is to ascertain the 
correct form of V peculiar to the particular case ; the form being found, the 
coefficients must be determined by experiment for each body : so far as theory 
is concerned, the mechanical classification of bodies would be complete, if the 
form of V were known for all. 

The conditions to be satisfied at the limits will be different, according as 
the surface is fixed, free, acted on by special forces, or in contact with other 
bodies. As there is no difficulty in forming the equations for any of these 
cases, I shall here give only the conditions at the limiting surface in contact 
with other bodies. 

Let F(x, y, z) = be the equation of the surface passing through the posi- 
tions of rest of the molecules which at any instant form the actual limiting 
surface /(ar, y, z, t) = 0, hence 

dF , dF J dF J A 

-j— dx + -j— dy + -T- dz = : 
dx dy * dz 

and if \, yu, v, be the angles which determine the position of the normal, we 
shall have 

dF 



cos \ = 


*~dx~' 


COS p = 


dF 
"dy 1 


COS v = 


dF 
dz 
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If w be an element of the surface 

a a dF 

ayaz = « cos A = «•« -=— , 

A A dF 

axaz = a cos il — kw -r— , 
dy' 

A A dF 

axay = w cos v = *•« -5— ; 
these equations will reduce the double integrals (4) to the form 

A = [[( d Z. d ^ + iZ. dF H d*-)^* 
JJ \da x dx da 2 dy da 3 dz ) 

rcfdV dF dV dF dV d£\ 
+ JJ \dfr ' dx + dfc ' dy + dfa ' dz J Ka ^ 

+ [[( d Z. d £+ d Z. d £ iZ.*H\ Kl ar 

JJ \dyi dx dy 2 dy dy s dz J ** 

These will be the double integrals resulting from one body, and from them 
must be subtracted similar terms derived from the body which bounds the one 
under consideration. F(x, y, z) — is common to both bodies when at rest, 
also £, ij, £, will be the same for all the bounding surface, f{x % y, z, t) = 0, during 
the motion ; also 8£, itj, 8£, are independent ; hence the condition at the limits 
will be 

A' - A" = 0, 

which is equivalent to three equations, 

/dV dV' '\dF rdV' Q dV' '\dF fdV' Q dVV\dF = 
\dai da x J dx \c?oj da 2 ) dy \da 3 da^ J dz 

dV' dV' '\d F (dV'o dV'jyF fdV' dV' '\dF 

rffr tfft J dx + \dp 2 dp 2 J dy + \dp 3 dfc ) dz ' *■ ' 

dV dV' '\dF fdV dV' '\dF /dV' dV' '\dF^ 
dyi dy! J dx \dy 2 dy 2 J dy \dy 3 dy s ) dz ' 

V', V" denoting the functions proper to the first and second body, and V 
denoting that the values of {x, y, z), deduced from F(x,y, z) = 0, have been 
substituted in V. 

To the three equations (6) must be added the self-evident geometrical 
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equations, which denote that the vibrating molecules at the bounding surface 
may be considered as belonging to either body ; they are three in number, 

& = 4o\ Vo = Vo\ & = &'• (7) 

Equations (6) and (7) contain the laws of reflexion and refraction of vibra- 
tions for all bodies, and are completely determinate when the form of Fis 
given for each of the bodies in contact. 

Equations (5), (6), (7), are necessary and sufficient to determine the propa- 
gation and reflexion of waves, so far as they are connected with each other ; 
and no mechanical theory of vibrations is correct which does not exhibit this 
connexion, or which assumes such laws of reflexion and refraction as contra- 
dict the laws of propagation ; the connexion between these laws is no proof of 
the truth of any theory, but the want of a connexion would be a proof of the 
inconsistency of a theory. 

The reduction of the general equations by the omission of the external 
forces may require some explanation. There are two kinds of waves, as I have 
stated in making the reduction, one only of which is the subject of our present 
inquiry. Fluid bodies, such as the atmosphere and the ocean, can propagate tidal 
waves depending on external forces ; they are also capable of propagating waves 
of sound which depend directly on the molecular forces ; solid bodies can only 
propagate the second species of waves. In considering this kind of vibration, 
we neglect the external forces, as they are of so much less intensity than the 
molecular forces, that they produce no effect on the motion ; but if we suppose 
the motion to cease, and inquire into the state of equilibrium of the body, we 
should then use the general formula, which includes external forces ; and, even 
in the case of motion, all the coefficients must be considered as variable, in con- 
sequence of the position of constrained equilibrium to which the body would 
return, if the motion ceased. 

It is evident from an inspection of equations (5) and (6) that the differen- 
tial coefficients of the function V, with respect to (a„ 02, &c), occupy an impor- 
tant position in the theory of elastic media. Hitherto, I have only given to them 
a mathematical definition, I now proceed to explain their physical meaning, by 
the consideration of an elementary parallelepiped of the body. If we conceive a 
plane drawn in any direction in the interior of the body, and consider the parts 
vol. xxn. p 
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of the body situated at opposite sides as acting upon any element of the plane ; 
in general, the effect of the particles at one side will be to produce a normal 
force, and tangential forces acting in the element ; these tangential forces may 
be resolved into two directions at right angles to each other. Let us now con- 
ceive an elementary parallelepiped, with one corner situated at the point (x,y,z) 
(these co-ordinates here denote the actual position of the molecule at any in- 
stant) ; let the forces acting on the side (dydz) be denoted by (P„ Q„ P,), P, 
being normal and parallel to the axis of (x), (Q,, P,) tangential and parallel to 
the axes of T and Z ; let the forces acting on the side (dxdz) be (P», Q 2 , P»)» 
Q 2 being the normal force ; and let the forces at the side (dxdy) be (P», Qs,B s ), 
P s being the normal force. The forces acting on the opposite sides of the 
parallelepiped will be : — on the side opposite to (dydz), 

( P, + -p dx J dydz, ( Q l + -— dx j dydz, ( P, + -j-± dx\ dydz ; 
on the side opposite to (dxdz), 

(p 2 +^dy\dxdz, (Q 2 + -^dy\dxdz, (R 2 + d ^dy\dxdz; 
on the side opposite to (dxdy), 

(Pj+ -si* dz) dxdy, (Q 3 + ~~dz)dxdy, (R 3 + -~dzjdxdy. 

These forces, acting on the six sides of the parallelepiped, must equilibrate the 
forces (Xdm, Ydm, Zdm), applied at the centre of the parallelepiped, and 
arising from external causes ; hence, the equations of equilibrium will be 






Zdm + [^P+^P + 



dR 2 dR*\ , , , « 
^ + ^dxdydz = 0; 

or, since pdxdydz = dm, and the molecular forces must act in the direction op- 
posite to the applied forces, including negative accelerating forces, 
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„ ,. dPi dPi dPt 
K X-u)=- d - + -^+ 1Fi 

fX r j\ dQ x dQa dQa /fl v 

r( Y -^=-dJ + -dJ+-d7'> (8) 

These are identical -with equations (2), and are true for all kinds of mole- 
cular action. In the particular case of a rigid parallelepiped, we should intro- 
duce another set of conditions arising from the equilibrium of couples. Let 
(Ldm, Mdm, Ndm) be external couples applied to the parallelepiped ; these 
must equilibrate the couples arising from the molecular action of the sur- 
rounding parts of the body ; it is easy to see that, neglecting the small forces 
arising from the differential coefficients of P lf P 2 , &c, the couples round the 
axes of ar, y, z, will be 

(-Ba + 0,) dydz, {P 3 + B,) dxdz, ( Q, + P,) dxdy. 

Hence the required conditions will be 

eL = i? 2 + Qj ; 
eM=P t + R l ; 
eN=Q l +P 2 ; 

and, if no external couples be applied, the conditions will be, 

P> = *, ; (9) 

Q l = P 1 ; 

since the couples (R it Q s ), &c, must act in opposite directions. 

These conditions (9) were given by M.Cauchy,* and afterwards adopted 
by M. PoissoN.f These writers seem to have considered them as necessary for 
all systems ; but this is not true, as equations (8) exhibit all the relations which 
exist between the forces and the motions produced. Equations (9) are necessary 

* Exercices de Mathematiques, torn. ii. p. 47. 

f Journal de l'Ecole Polytechnique, cahier xx. p. 84. 

p2 
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for the equilibrium of a rigid parallelepiped, and will be shown in this memoir to 
be satisfied by the equations of equilibrium of bodies whose molecules attract and 
repel each other in the direction of the line joining them ; but if no supposition be 
made as to the nature of the molecular action,there will be no condition resulting 
from the equilibrium of couples ; for we have no right to assume, in the equili- 
brium of a parallelepiped, whose elements may alter their relative position, and 
thus develope new forces, that the same equations hold as in the equilibrium of 
an isolated rigid parallelepiped, for which case six equations of condition are 
necessary, arising from the equilibrium of forces and of couples. 

Restoring the usual signification of (x, y, z) in equations (8), and omitting 
the external forces, we obtain 

d^_dP x dP 2 dP t 

d(? ~ dx dy dz ' 

__ d*r,_dQ } dCL dQi 

dt* ~ dx* dy + dz ' V V) 

_ e <££_dR 1 dRj dRt 

dt* dx dy dz ' 

These equations correspond with (5), and by comparing them we may deduce 
the following relations : 

n= P% d z=p 3 , d z=p. 

da x ' ' da 2 ' da 3 3 ' 

dV_ n dV_ n dV_ n 
dfr-^' dp,'*' dj 3 ~^' 
dV dV_ p dV 

dy, dy 2 dy 3 

from which we obtain the following theorem : 

" If through any point (a?, y, z), three elements of planes be drawn parallel 
to the co-ordinate planes, the total action of the part of the body lying at one 
side of these planes will consist of three normal and six tangential forces ; and 
these forces may be expressed by the differential coefficients of the function V, 
with respect to (a,, aj, 03, &c.)." This theorem is of importance in classifying 
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bodies, as it enables us to pass directly from the resultant forces of the mole- 
cules to the form of the function which determines the laws of propagation 
and reflexion. Introducing these forces into the conditions at the limits (6), 
we obtain 

p,dF dF dF_ dF d_F dF 

dF dF dF_ dF dF dF^ 
^fa+^fy + ^fo-^te+^dy+^dz' 

dF dF dF _ dF dF dF 
jK "^ + m dy~ + iqa dz~^ l d^ + °*4y °*dz~' 

If the axis of z be made to coincide with the normal, we shall have 

^ = ^ = ^=1. 

dx ' dy ' dz 

Hence, 

Pi _ jyn w %ji 

OS — * 03) So — so > 

<&=<&, rt = ijl\ (li) 

Jfl Till yl _ yll 

*H% — -"oa i &o — &o • 

These equations at the limits are the mathematical statement of two facts, of 
which one is mechanical and the other geometrical. 

1. That the forces, normal and tangential (arising from molecular action), 
acting upon an element of the bounding surface, must be equal and opposite for 
the two bodies in contact. 

2. That the motion of the particles in the bounding surface may be con- 
sidered as common to both bodies. 

I shall now return to the general equation (4) ; let the function Fbe di- 
vided into homogeneous parts, so that 

V- 00 + tf>, + 02 + & + &C-. 

where </> is constant, 0, of the first degree with respect to (oj, a 2 ), and so of the 
others; then, neglecting <f> and terms higher than the second, we obtain 
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0, = Aft + A 2 02 + Aoj 

+ B& + B 2 p 2 + B& 
+ C 1 y l +C 2 y 2 +C i y i . 

fc = (a 1 s )a 1 ' + (a 2 *)a 3 J +&c. 
+ {a x bi) • a x p 2 + &c, + (c 3 a 2 ) • 7,0, + &c. ; 

the expressions within the brackets denoting the corresponding coefficients. 

fa will contain nine coefficients, which will be constant, or functions of 
(x,y,z), according as there are no external forces, or vice versd; and <fc will 
contain forty-five distinct coefficients. I shall consider these functions separately. 
Introducing fa for Fin equation (4), we obtain, 

III 6 d? l ^W 8l7+ 3 8f ) dxdyd * = M A &+ S fo+ <W *y**> 

+ \\(AM + B 2 lri + C£g) dxdz, ( 12) 
+ IJ(AM + Bfr, + Cft) dxdy 



whence we obtain 



d 2 £ __ <L4 t «L4 2 d4 3 
df~dx dy dz ' 



<f ^ ^A dLB, d2? 3 / . 

£ ^ = d^, dC.dC, 
d? ~~ dx dy dz ' 

The terms on the right hand side of equation (13) must be added to the 
dynamical equations arising from fa, even in cases where (X, Y, Z) may be 
neglected on account of the intensity of the molecular forces. If we wish to 
take account of the external forces, we should use, instead of (13), the fol- 
lowing : 
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d*i- _dA, dAj dA 3 „ 
d?~ dx dy dz ' 

C df ~ dx + dy + dz ; 

<P£_dC } dC t dC 3 „ 
d£~ dx dy dz 

If no external forces act, equations (13) disappear, in consequence of A ly A it 
&c., becoming constants; but the conditions at the limits (12) will still remain, 
unless the coefficients be not only constant, but zero. 

The part of the differential equations of motion depending on fa will be 
found by introducing 8^ 2 in place of iV in equation (3), and integrating by 
parts, according to the methods of the calculus of variations ; but it may be 
found more readily by introducing fa for V, in equation (4). Neglecting the 
equations of condition at the limits, we obtain for the equations of motion, 

~ e dJ = ^d? + ^df + ^d? 



<*V , . N A , , ..d'ri 



+ («*)S+(-*)^l+(«A)5J 

+ (aA + a A) ^ + (aA + a A) gjj- + (<*A + «A) ^ 

+ («iC)g+(aA)0 + (^)g 
+ (a 2 c J+ a sC2 )^+(a lCs +a3C 1 )^+(a 1 c 2 + a 2 c 1 )^. 






(14) 
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+(»a)3+(w§+(»*)S 

v72j* J%t* r&t* 

+ (^i) a? +(^)^+(fA) a? 

A%£ J2t* Alt' 

+ (a* + «*) ^+ (oa + oa) jgj + (a.4 + a*) j^ 

By combining these equations with (12), and comparing (4), we obtain from the 
function V = ^ + #2 the following equations of motion, 

ar dxydc^ J dy\da2 / dz\da 3 ) 

-•§+=($+*) + *(%♦*)♦=(£♦*) (15) 
_.g <(* + c,Vf(£+ftV£(*+<0 

These are the equations of small oscillations of all media, for which the direct 
influence of the external forces may be neglected. 

In equations (14) it will be observed that nine of the coefficients are 
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composed of two terms, coefficients of the original function 2 > these compound 
coefficients will, in the equations of motion, be equivalent to simple coefficients, 
so that the total number of distinct constants in the equations of wave-motion 
will be reduced to thirty-six ; this reduction cannot, however, be introduced 
into the conditions at the limits, because each coefficient will at the limits be 
multiplied by a different function of (x, y, z). Hence may be deduced an im- 
portant theorem, which I shall prove in a general manner. 

" Two bodies, different in their molecular structure, may have the same 
laws of wave-propagation, but cannot have the same laws of reflexion and 
refraction." 

The laws of wave-propagation (5) depend upon the differential coefficients 

of f -j— , &c. j, taken with respect to (x, y, z) ; but the laws of reflexion and 

refraction (6) depend on the quantities f -j-, &c. } themselves. Each of these 
quantities is of the form 

dx dy dz dx dy dz dx dy dz 

If, therefore, for example, in the first of the equations (5), there be in (-7- ) 

a term of the form LB -A there will be in -r- a term of the form \A'-$ 
\ dy) da* \ dx 4 

then, differentiating the first with respect to x, and the second with respect to 

d*£ 
y, we shall obtain, as part of the equations of motion, (B + A') ~- , while 

cbxcLy 

the corresponding part of the conditions at the limits will be [[ B -~ l^dydz 
+ \\A' -j-fydxdz. If, now, we suppose another body such that (A'~) and 



dy) 

[B-y-J are terms in (-j—) and (-r-J, the part resulting from these quan- 

d 2 £ 
tities in the equations of motion will be (A' + B) -~- , which is identical with 

(LX&'U 

the former value, while the corresponding portion of the conditions at the 

VOL. XXII. Q 
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limits will be I -4' -3- $%dydz +\\B-j- l^dxdz ; which cannot be reduced to the 

former value, unless A' — B. 

It is sufficient here to establish the general theorem. I shall return to it 
again, in comparing Professor Mac Cullagh's theory of light and that ad- 
vocated by M. Caucht and Mr. Green, which give the same laws of wave- 
propagation, but differ in the laws of reflexion and refraction. 

I shall now integrate the equations of motion (14), for the particular case 
of plane waves and rectilinear vibrations. Let (/, m, n) be the direction cosines 
of the wave-normal, (a, /3, 7) the direction of molecular vibration, and (») the 
velocity of the wave ; then the integral will be 

§ = cos a *f\ —(lx + my + nz — vt) >, 

r\ = cos /3 •/ -j — (h + my + nz — vt) >, 

£ = cos 7 •/ 1 -£ (1$ + my + nz- vt) \. 

Introducing these values into the equations (14), we obtain 

ev 2 cos a = F cos a + £F cos /3 + Q' cos y, 

ev 2 cos j3 = Q cos j8 + F' cos 7 + H' cos a, (16) 

ev 2 cos 7 = R cos y + G' cos a + F cos /3 ; 
where 

P' = (a, 2 ) I 2 + (a/) m 2 + (a 3 2 ) n 2 + 2(a 2 a 3 ) rnn + 2(«ia 3 ) In + 2(a l a 2 ) Im ; 
Q' = (&, 2 ) I 2 + (h 2 ) m 2 + (b 3 2 ) n 2 + 2(&A) m + 2(&A) fti + 2(* 1 * 1 ) An ; 
^' = (Ci 2 ) P + (c 2 2 ) to 2 + (c s 2 ) ri* + 2(e 2 c 3 ) w» + 2(c,c 3 ) fo + 2(c,<? 2 ) 2m; 
F' = (^iC,)/ 2 +(t2C2)m 2 +(i 3 c a )n 2 +(i 2 c s +63C2)mn-i-(5iC 3 +63Ci)^n+(V2+&2Ci)/m ; 
£' = (aiCi)P+(a 2 c 2 )m 2 +(a 3 c 3 )n 2 +(a 2 c 3 +a 3 c 2 )mn+(a x c 3 +a 3 c x )ln+(a i c 2 +a 2 c l )lm ; 
R'zs. (a,6, ) l 2 + (a A) m 2 + (a 3 6 3 ) n 2 + (a 2 6 3 + 03^2) wm+ (a^a+aa&i ) &»+ («i #2+ a A ) fon. 

Equations (16) are the well-known equations for determining the axes of 
the ellipsoid whose equation is 

Fx 2 + Q'y 2 + Ez 2 + 2Fyz + 2 0'aw + 2£T^ = 1. (17) 
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There are, therefore, three possible directions of molecular vibration for a 
given direction of wave plane ; and there will be three parallel waves moving 
with velocities determined by the magnitude of the axes of the ellipsoid} the 
direction of vibration in each wave being parallel to one of the axes. The 
theorem involved in equation (16) was first proved by M. Cauciiy * for bodies 
whose molecules act by attractions and repulsions in the line joining them. It 
is here extended to every kind of molecular action, and shown to be a funda- 
mental property in molecular dynamics. It may be worth while to examine 
the reason of its truth. It arises from the homogeneity of equations (14) 
resulting from the absence of external forces. The right hand member of each of 
the three equations of motion consists of eighteen terms; there will, therefore, 
be in all fifty-four terms ; if each of these were supposed to have different and 
independent coefficients, equations (16) would cease to be true, and the coeffi- 
cients of (cos o, cos j3, cos 7) would be nine distinct quantities, so that the 
theorem which represents the direction and magnitude of the molecular vibra- 
tion by means of an ellipsoid whose coefficients are functions of the direction 
of the wave, would no longer be applicable. Equations (14) contain only 
thirty-six distinct constants, and this reduction in the number of constants arises 
from the assumption that the virtual moments of the system may be repre- 
sented by \\\l Vdxdydz ; which is equivalent to assuming that the virtual mo- 
ments of the molecular forces applied at any point may be represented by the 
variation of a single function : 

Qlq + Q'lq 1 + &c. = I V. 

The cubic equation whose roots are the squares of the reciprocals of the 
axes of the ellipsoid (17) is 

(F-s)(Q'- s )(R- s )-F' i (F-s)-G'*{Q'- s )-IT*{R f - s )+2F'G'H , =:0; 

where s=ev\ Hence, if (P, Q, R, F, G, H) denote the same functions of (x, y, z) 
that {P, Q', R, F', G', £F) are of (I, m,n), it maybe shown, in a manner similar 

* Exercices de Mathematiques, torn. v. p. 32. 

Q2 
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to the method for bodies whose molecules act in the line joining them, that the 
surface of wave-slowness* will be 

(P-l)(Q-l)(B-l)-F*(P-l)-G 2 (Q-l)-H 2 (R-l) + 2FGH=0. (18) 

(P— 1, Q — 1, &c.) equated to zero, denoting the six fixed ellipsoids used in 
the memoir referred to ; these ellipsoids thus appear in the general theory of 
molecular dynamics, and will have the same use in- interpreting the conditions 
at the limits in the general problem, as in the particular case of attracting and 
repelling molecules. As their use has been fully explained in my former me- 
moir, I shall only refer to it, and pass on to other subjects. 

I shall here state the general conditions necessary, in order that the mole- 
cular equilibrium of a body removed from the influence of external forces 
should be stable ; this investigation will be found of use in the subsequent part 
of this memoir. The equation of virtual velocities, in the case supposed, will 
become 

or, as it may be more accurately stated, 

\\\ZVdxdydz = 0, or < ; 
this sum never becoming positive for possible displacements. The equa- 
tion of virtual velocities, as stated by M. Poisson and other writers, supposes 
no virtual displacements but those for which equal and opposite displace- 
ments are possible ; the correct statement of the principle is, that the sum of 
the virtual moments of the system can never become positive for possible dis- 
placements. For the full development of this important correction of the equation 
of virtual velocities, as given by Lagrange, I shall refer to a memoir of M. 
Ostrogradsky, contained in the Memoires de VAcademie de St. Petersbourg.f 

The application of the principle to the present case will be evident upon 
the statement of the question. As the form of Fis given in terms of (a, , 02 , &c), 
we shall have 

by — -j— ba, + -=— 8o 2 + &c. ; 
aai ddi 

* Vide Transactions of the Koyal Irish Academy, vol. xxi. part ii. p. 172. 

f Tom. iii. p. 130. 
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fdV dV \ 

(-j— , -j— , &c. ) denoting forces tending to alter the quantities (a,,a 2 ,a 3 ,&c.), 

\dfiY da 2 J 

which are functions of (x, y t z). The whole body may, therefore, be divided 

into couches by a series of surfaces whose equation will be 

a,-C = 0; 

C denoting the parameter of the system. In a similar manner, the body may be 
conceived as divided into couches by other sets of surfaces corresponding to 
(a 2) a 3 , &c). In any one of these couches the corresponding function 
(a u an &c.) will have a constant value, which will vary from one couche to 
another. 

If L = be the equation of a surface, we may conceive all space as divided 
into two portions, which will be distinguished by the property, that the portion 
lying at one side of the surface will have the function of (x, y, z) denoted by 
(L), positive ; while for the rest of space, lying at the other side of the surface, 
the function (_L) will be negative. Similarly, 8i will be positive for all dis- 
placements made at one side of the surface ; negative for all displacements on 
the opposite side ; and zero for displacements in the surface itself. 

Let us now resume the equation, a, — C = 0, which denotes a surface drawn 

rlt 
in the interior of the body, along which the value of ~ remains constant. It is 

necessary for the stable equilibrium of the body, that if the particles composing 
this surface be displaced from it, the molecular forces developed by the dis- 
placement should tend to restore them to the surface ; this condition will 

(dV dV \ 
require that I -7— , j— , &c. 1, which are the forces developed by the displace- 
ments (8a,, 802, &c), should have signs opposite to those of the displacements. 
If no forces are developed tending to restore the molecules, we shall have 

dV_ a dV 
So"" ' Sa7 = °' &C - 

and these equations will determine the limits of stable and unstable equilibrium 
in the body. 
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SECTION II. — LAWS OF NOKMAL VIBBATIONS. 

It has been shown that the differential coefficients of the function V, taken 
with respect to (ctj, dj, &c), denote the normal and tangential actions of the 
surrounding body on an elementary parallelepiped. I shall now suppose that 
the body is so constituted that the tangential forces vanish, and so deduce the 
laws of a body capable of transmitting pressure in a normal direction only. 
The condition that the tangential forces vanish will give the equations, 

<LZ = o <*I-o ^-o 

da 2 V ' 48, ' dy t ~ "' 

d -Z=o d ~l=o d -l=o 

da 3 dp 3 ' dy 2 ' 

which will reduce the function to the form 

2 V = A a, 2 + Btf + Cy* + 2Pj3*y 3 + 2Qo,7 3 + 2Ra& . 

At first sight, it might appear that we might assume this or any other form 
of function to define a body, and proceed to deduce the laws of motion ; but, 
as it is possible that an assumed form of the function may be only true for par- 
ticular axes of co-ordinates, it is always necessary to examine whether a trans- 
formation of co-ordinates will introduce any new terms ; if this be the case, 
then the assumed function will not represent completely the molecular structure 
of the body, as it will be merely a simplification of a more complex function, 
produced by assuming particular axes of co-ordinates, which are connected 
with the crystalline structure of the body. I shall examine the function just 
given for normal pressures by this method, and determine it so as to be inde- 
pendent of the axes of co-ordinates. The relations between the two systems 

of variables are 

x = ax' + byf + cz 1 , 

y = a'x' + b'y' + c'z', 

z = a"x'+b"y'+c"z'; 

also (£, 77, £) will satisfy these equations. Let us now assume 
« = ft + 72i w=7i + «3» w = o, + 0„ 
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we may immediately deduce, by the formulas for transformation of the indepen- 
dent variables, the following values for the nine quantities (o„ 02, &c.), 



a, = aV, + & 2 $ 4- <?y 3 + bcu' + acv' + abw', 

fr = a' 2 a\ + b'% + c' 2 y' 3 + b'c'u' + a'c'v' + a'b'w', 

y 3 = a" 2 a[ + b"% + c" 2 y' a + b"c"u' + a"c"t/ + a"b"w' 



r a' (a"a\ + b"^ + c"^) 
fr= \ + b'{a"p[ + b"fr + c'%) 
+ c'(a"y' l + b"y' i + c"y 3 )\ 

■ a" (aa| +ba i + ca 3 ) 
+ b"(ap\ + bfc + cp 3 ) 
+ c"(ay' l + by' i +cy 3 ). 

r a (a'a[ + i'aj + cX) 

a,= \+b(a'^ + b'a^ + e'^' 3 ) 

[ + e(a'y' l + b'y' i + c'y' 3 )_ 



7i = 



■ a" (a'a[ + b'<£ + c'cQ] 
7 ,= + 6" («'# + &'# + <;'&) 
+ c"(aV. + ^ + ^)j 

r a (o'V + V'cd, + c"a 3 y 
a 3 =\+b(a"ft + b"fi' 2 + c"fi' 3 ) 

l+c(a"y[ + b"y' a + c"y' 3 )_ 

(" a' (aa'i + ba^ + Ca 3 ) 

ft =+*>'(«#+&# + <$) 
[ + c'(ay\ + by' i + cy 3 ) 



(19) 



It is plain from the first three equations that the change of the co-ordinate 
axes will introduce into the function V three new variables (u, v, w). Hence, 
the function which I have assumed is only admissible for particular axes, and 
the true form of the function from which it is derived will be that of a function 
of the six variables used in my former memoir. 

It is, however, possible to obtain such a function of ( -5-, -^, -£ ) as shall 

not change its form with the transformation of co-ordinates. For, adding toge- 
ther the values of these quantities, we shall banish («', v', uf), on account of the 
relations, 

be + b'c' + b"c" = 0, 

ac + a'c' + a"c" = 0, 
ab + a'b' + a"b" = 0. 

The result of the addition will be 

dx dy dz dx' dy' dz' 
This quantity («) will retain its value, independent of the directions of the 
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co-ordinates, and consequently the proper form of the function V, for bodies 
which transmit only normal pressure in every direction, will be 

V= F(w). (20) 

There are two other forms of the function V (considered as containing the 
nine differential coefficients) which possess this property of reproducing them- 
selves by transformation of co-ordinates. Let (X, Y, Z) be determined by the 
following equations, 

X = /3s — 72, Y=y t — a 3 , Z — a^ — ^x. 

It appears immediately from equations (19) that (X, Y, Z) are expressed by the 
following relations in terms of (X', Y', Z'\ 

X=aX'+bT+cZ\ 
Y=a'X'+b'r+c'Z', 
Z = a"X' + b"Y' + c"Z'. 

These formulae are well known, and prove that another self-producing form of 

the function V will be 

V=F(X, Y,Z). (21) 

I have proved in my former memoir that if (a, 0, 7, u, v, w) denote the co- 

efficientS {§' Tf % (^ + |) (I 4 "!) (f + S) W™ UbeCOme « 
by transformation of co-ordinates, linear functions of (a, ft, y, ut, v', w') ; and 
that the formulae of tranformation are the 6ame as for (tf 2 ), (j/ 2 ), (z 2 ), (2yz), 
(2xz), (2xy). Hence, a function of these six quantities will reproduce itself, 
by transformation of axes of co-ordinates. The third form of V which possesses 
this property is, therefore, 

V=F(a,p,y,u,v,w). (22) 

The forms of the function V, which have just been determined, are perfectly 
general, and do not merely express properties of the body with reference to 
particular axes, but general properties of its molecular structure, independent 
of the directions of the co-ordinates. 

I shall consider in the present section the function (20),which is peculiar 
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to bodies which transmit normal pressure in every direction. Before deducing 
the equations of motion, it may be useful to show, in another manner, the neces- 
sity for reducing the function of (a, /3, 7) to a function of (a + j3 + 7). By the 
theorem proved in my former memoir, (a, /3, 7) are transformed by the same 
equations as (x 2 , y*, z 2 ) ; hence, V= F(a, /3, 7) may be represented by a surface 
whose equation contains only the squares of the co-ordinates. It is evident 
that such a surface will be symmetrical with respect to the co-ordinate planes ; 
and if this condition be satisfied for every system of co-ordinates, the surface 
must be a sphere ; hence, 

V = F(a,(3 1 y) = F(a + f} + y ). 

The equations of motion deduced from (20) will be the equations commonly 
used in hydrodynamics, and for this reason it may be useful to state them in 
their most general form. Let (x, y, z) denote, as in equations (1) and (2), the 
actual positions of the molecules ; then, since 

\dx dy dz J 
we shall obtain from equations (2) 

1 dp „ du du du du 

p ax at ax dy dz 

1 dp -v dv dv dv dv /na . 

- -f-= * --r-u j v -j w-r-; (23, a) 

p dy dt dx dy dz 

1 dp r? dw dw dw dw 

- -f- = <6 — -5- - u -j u-5 w-j- ; 

p dz at dx dy dz 

dV 
assuming p = -j-, and recollecting that (u f , v', w') are the total differential co- 
efficients of (u, v, w), taken with respect to (t). These are the mechanical 
equations of the problem, and, combined with the equation of continuity, 

ffe+^ff) + 1M + ^£l = 0, (24, a) 

dt dx dy dz 

vol. xxn. b 
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contain all that is requisite to determine the motion of points situated in the 
interior of the mass. The mechanical conditions to be satisfied at the limits 
are contained in the double integrals, which become, by the substitution of the 
element of the bounding surface, 

l\pdS (Zx cos \ + Zy cos ft + Zz cos m), 

(A, fi, v) denoting the direction of the normal. It may be easily shown from 
this expression, that if P denote a function of (x, y, z, i), expressing the normal 
forces applied at each point of the surface, the mechanical condition at the limits 
will be expressed by the equation, 

p-P = Q. (23, b) 

To which must be added the equation of the bounding surface, deduced from 
geometrical considerations, 

dF dF dF dF A /0 , . v 

-rr + u-r- + v-j- +w-i-= 0. (24, b) 

at ax ay dz 

It is unnecessary to enter further into this subject, as it is fully treated by 
Lagrange, in the second volume of the Mecanique Analytique, and is indeed 
the only example given by him of the application of his formulae to the problem 
of bodies composed of continuous points. 

Resuming the former signification of (x, y, z, £, r\, £ ), we find from equations 

dV 
(5) or (15), since 2V = 2pw + Aa> 2 , and -j-=p + Aw, 

d 2 £ _ dp . a ^ 
dt 2 ~ dx dx' 

cPq _ dp . dw 
~ e dF~dy + dy' 

d 2 g _dp . dw 
d? ~ dz dz' 

If no external forces act, p and A will be constants ; and, in this case, these 
equations become, assuming a negative sign for A, so that the equilibrium may 
be stable, 
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d*% A dt* 



*d?~ A dx' 



d 2 v _ A dm (25) 



*d?~ A dz 



and, by differentiating with respect to (#, y, z), we obtain 

d i m __ A fd?i» d?a> d 2 w\ ,„e\ 

rf?~7^ + ^ 3+ '3iy , { } 

which determines the cubical compression as a function of (#, y, z, t). 

The equations (25) and (26) are applicable to the propagation of sound in 
gases, and to the longitudinal vibrations of elastic solids ; but there will be an 
essential difference between the two cases. If the gas were freed from the 
action of external forces, and unconfined at its limits, its molecules would sepa- 
rate, and a displacement would develope no force tending to restore them to 

their original positions ; hence, l-j- iw J would become positive, which is in- 
consistent with stable equilibrium, and the velocity of wave propagation 
( nj—} would become imaginary. In order, therefore, that the equilibrium of 

the gas be stable, we must suppose a constant pressure exerted at the limits, 
sufficient to keep the molecules together, and restore them to their original 
positions, if displaced. No such condition is requisite in the elastic solid, for 

i-T- lw j will be negative, without the assistance of forces applied at the bound- 
ing surface. 

SECTION III. LAWS OF TBANSVEESB VIBEATIONS. 



In order to obtain the form of the function V peculiar to transverse vibra- 
tions, we must suppose that, if a plane be drawn through the body in any 
direction, the molecular forces exerted on any element of the plane will be 
altogether tangential. Hence we obtain 

r2 
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dV . dV n dV 
aaj aft dy 3 

and the function V will become 

F=.F(ai, a,, ft, ft, 7k, 72)- 

But, as I have already shown, such a form must be assumed, as will reproduce 
itself by transformation of co-ordinates. The only function which satisfies this 
condition is the function (21), 

V=F(X,Y,Z). 

This function, deduced in a different manner, has been used by Professor 
Mac Cullagh in his mechanical theory of light ; and for the discussion of its 
properties and the laws of propagation, reflection, and refraction, deduced from 
it, I shall refer to his memoir in the Transactions of the Royal Irish Academy.* 
As, however, I shall have occasion to use them hereafter, I shall here state the 
differential equations of motion, and the conditions at the limits. On account 
of the form of the function, the following relations exist : 

dV = dV = _dV 
dX dft dy 3 ' 

dV = dV ;= _dV 
dY dyi da z ' 

AV = dV = _dV 

dZ da.2 dft } 

Hence equations (5) will become 

d*Z_d_ dV_£ dV 
e df~dy'dZ dz' dY' 

d\ d_ dV_d_ dV 

dz ' dX dx' dZ' ™"> 



e 'df 



d%_d L dV_d L dV 
S di*~dx'dY dy'dX' 



Vol. xxi. p. 17. 
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These are the equations of motion; also equations (6) will become 



dV' a dV\ 







' '\dF fdV dV' '\dF_ 
'J dy \dY dY ) dz ' ' ' 



dZ dZ 

(dV dV' '\dF (dn_dVT\dF_ 
\dX~~dXjdz~~\dZ dZjdx-"' ^° } 

fdV _ dV' ' \ dF _ fdj^ dV'o'\ dF_ 
\dY dY) dx \dX dXj dy 

These are the conditions to be fulfilled at the limiting surface F(x,y,z) = 0. 
They are equivalent to two conditions only, as the third equation may be de- 
duced from the first two. The reason of this is evident a priori; for the con- 
ditions at the limits express in general, that the normal and two tangential 
pressures arising from the molecular forces in each body equilibrate each other 
for every point of the separating surface ; but in the present case there are no 
normal pressures ; hence there can be only two mechanical conditions at the 
limits. To these conditions at the limits, arising from the mechanical equa- 
tions, should be added the three geometrical conditions resulting from the equi- 
valence of vibrations ; these conditions are 

So = So > Vo = Vo j £ o = i o • 

These, together with the mechanical conditions, will be equivalent to five equa- 
tions, which Professor Mac Cullagh reduces to four by the hypothesis that the 
density of the luminiferous medium is the same in all transparent bodies ; this 
hypothesis is necessary in order to reduce the number of equations to the num- 
ber of unknown quantities in the problem. 

If no external forces act upon the system, the function V will be reduced 
to a homogeneous function of the second order, 

2 V = PX 2 + Q Y 2 + RZ 2 + 2FYZ + 2GXZ + 2HXY; 

and since (X, Y, Z) are transformed by a change of co-ordinates, in the same 
manner as {x, y, z), it is evident that the coefficients of the rectangles will 
vanish for axes of co-ordinates which coincide with the axes of the ellipsoid, 

Px 2 + Qy 2 + Rz 2 + 2Fyz + 2Gxz + 2Hxy = 1. 
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Hence the simplest form of the function will be 

2V=PX*+QY* + RZ\ 
The axes of co-ordinates are in this formula the axes of elasticity. 



SECTION IV. EQUATIONS OF A SYSTEM WHOSE MOLECULES ATTRACT AND REPEL 

EACH OTHER. 

The third form of the function V, which possesses the property of repro- 
ducing itself by a transformation of co-ordinates, is given by equation (22), 

V=F(a, /3, 7, u, v, w). 

The six variables contained in this function are the quantities upon which a 
change in the distance between the molecules depends ; hence the variation of 
this function will express the virtual moments of a system of forces tending to 
alter this distance. In my former communication* I have made use of this func- 
tion, using definite integrals to represent the coefficients ; the equations thus 
found contain a smaller number of constants, than if the coefficients had been 
assumed to be' arbitrary, without any relation expressed by definite integrals. I 
shall here use the function in its most general form, as the use of definite inte- 
grals may, perhaps, involve an hypothesis which would be too restricted to 
represent all the bodies whose molecular forces act in the line joining the 
molecules. 

The following relations will exist, in consequence of the form of the function, 

dV = dV = dV 

du ~ dfo ~ dy 2 ' 

dv dyi d<h ' v ' 

dV_dV_dV 

dw ~ rfa 2 — d$ x ' 

These equations will reduce the equations of motion (5) to the following : 

* Transactions of the Royal Irish Academy, vol. xxi. p. 151. 
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(30) 



These are equivalent to equations (24) of my former paper, and are an imme- 
diate consequence of (29), which express the relations among the resultants of 
the molecular forces consequent upon the restricted form of the function. Equa- 
tions (29), or the corresponding equations (9), 

0,-^ = 0, R 1 -P 3 = 0, P 2 -Q, = 0, 

are analogous to a statical theorem given by Fbeskel in his Memoir on double 
Refraction. If no external forces act upon the system, or if their influence may 
be assumed to be constant at all points of the body, then the coefficients of V 
will become constants, and the linear part of V will introduce no terms into 
the differential equations, which will in this case depend upon a homogeneous 
function of the second order. There will be, therefore, in general, twenty-one 
constants, viz., the coefficients of the function 

2V= (a> 2 + (/3 2 )/3 2 + (y 1 )^ + (u 2 )u 2 + (»>* + (w 2 )w 2 

+ 2{ (^7)^7 + (<»Y)a7 + (a/3)a/3j + 2\(vw)vw + (uw)uw + (uv)uv\ 
+ 2uj(au)a + (/3tt)/3 + (yu)y\ + 2v{(av)a + (/3t>)/3 + (7^)7} 
+ 2wj (aw)a + (/3w?)/3 + (7^)7}. 

The equations of motion deduced from this form of the function V are 
identical with the equations used by M. C auchy in his theory of light.* They 
are deduced by M. C auchy directly from the consideration of attractive and re- 
pulsive forces between the molecules. These equations have been used also 
by Mr. Green, who seems to have considered them as more general than M. 
Catjchy's equations.! 

The function (22) is not altered in form by an alteration of the co-ordinate 

* Exercices de Mathematiques, torn, v., p. 19. 

t Cambridge Philosophical Society's Transactions, voLviL p. 121. 
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axes, and is the most general which can be used for the case considered in the 
present section. The direction of the molecular vibration will not be either 
normal or transversal, but will be determined by the axes of an ellipsoid first 
noticed by M. Catjchy in treating of bodies whose molecules attract and repel 
each other. It may be worth while to examine whether there are any subordi- 
nate functions which possess the property of not being changed in form by the 
change of co-ordinates ; as such functions, if they exist, will contain the laws 
of classes of bodies contained under the general function (22). I have already 
discussed one such function (20), which is evidently a particular case of (22). 
I shall now determine another remarkable subordinate function, which, while it 
is much less general than (22), yet contains most of the bodies whose proper- 
ties are at all interesting. The following formulae of transformation may be 
easily demonstrated.* 

a = a 2 a! + b 2 fi + &y' + bcu' + acv' + abw'; 

p = a'V + b' 2 p' + c'V + b'c'v! + a'c'v' + a'b'w'; 

7 = a" 2 a' + b"Y + c"V + b"c"u' + a"c"v' + a"b"w'; 

u^2aW'a'+2bV'^+2cY'y , +(b'c , '+b'Y)u , +(aY'+a''cy'H^''+a%')w'-, 

v = 2aa"a'+ 2bb"p' + 2cc"y'+ (bc"+b"c)u'+ (ac"+ a"c)v'+ (ab"+a"b)w'; 

w = 2aa'a! + 2bb'ft + 2cc'y + (be' + b'c)u' + (ad + a'c)v' + (ab' + a'b) w'. 

If we now assume six functions of (a, /3, y, u, v, w), defined by the following 

equations : 

\ = u 2 — 4/3y, fi = v 2 — 4ay, v = W 2 — 4a/3 ; 

= 2au — vw, x = 2j8v — uw, ■$■ = 2yw — uv ; 

it may be shown without much difficulty, that these new functions are trans- 
formed by the following equations : 

\ = a 2 \' + &V + cV + 2bc</>' + 2ac x ' + 2ab^'; 

fi = a' 2 \' + 6'V + c'V+ 26V*' + 2aV X ' + 2a'b'+'; 

v = a" 2 X' + 6'V + c"V + 2b" c"^ + 2a"c" x ' + 2a"b"^' ; (31) 

= a / a"\'+b'b"n'+c'c"v , +(b'c"+b"c')<l>'+(a'c"+ a"c')x'+( a ' h "+ a " b ')^'^ 

X = aa"\' + bb"(x' + cc"v' + (be" + b"c)<f>' + (ac" + a"c)x' + (ab" + a"b) f; 

f = aa'X' + bb'fi' + cc'v' + (be' + b'c) <f>' + (ac' + a'c)x' + (ab' + a'b) +'. 

* Transactions of the Royal Irish Academy, vol. xxi. p. 160. 
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Hence a function of (X, p, v, </>, x> "^) will reproduce itself by transformation 
of co-ordinates ; let the function be 

27= P\ + Qfi + Rv + 2F(f> + 2G X + 2#f 

I shall first prove the existence of three rectangular axes, for which this 
function reduces to its first three terras. If the axes of co-ordinates be trans- 
formed, and the coefficients of (<ft , %', ^') equated to zero, we obtain 

Pbc + Qb'c' + Rb"c" + F(b'c" + b"c') + G(bc" + b"c) + H{bc' + b'c) = ; 
Pac + Qa'c' + Ra"c" + F(a'c" + a"c') + G(ac" + a"c) + H(ac' + a'c) = ; 
Pab + Qa'V + Rx"b" + F{a'b" + a"b') + G(ab" + a"b) + H(aV + a'b) = 0. 

These equations will be satisfied by assuming for axes of co-ordinates the axes 
of the ellipsoid 

iV + Qy 2 + Rz 2 + 2Fyz + 2Gxz + 2Hxy = 1. 
Hence, for these particular axes, 

2 V = P(u 2 - 4/*y) + Q(f - 4a 7 ) + R(w 2 - 4a/3). (32) 

Using this value of Fin equations (30), we obtain for the equations of motion, 

d^_ dZ dY 

e de- K d y ~ u dT> 

d i t ] _ p dX dZ ,„„. 

_^_ Q dY_ p dX 

dt 2 dx dy 

(X, Y, Z) denoting the same functions as in (27). If the function were 

2V=PX 2 +QY* + RZ 2 , (34) 

equations (27) would be the same as (33). 

These equations are those used by Professor Mac Cullagh and by Mr. Green. 
They denote transverse vibrations, and will give Fresnel's wave-surface for 
plane waves, and also the vibrations parallel to the plane of polarization. Mr. 
Green has deduced his equations from a homogeneous function of the second 
order of (o, /3, y, u, v, w), by restricting it so as to be capable of propagating only 
VOL. xxii. s 
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normal and transverse vibrations; this restriction will reduce the constants 
from twenty-one to seven; six of which belong to the transverse vibrations, 
and the seventh to the normal vibration. It is to be remarked, however, that 
a function of (\, /*, v, 0, %, ■$■) will represent only a part of the properties of 
bodies whose molecules attract and repel each other. Such a body is always 
capable of transmitting normal vibrations, and though it may transmit trans- 
verse vibrations following the laws of Fbesnel's wave-surface, yet the normal 
vibration cannot be supposed to vanish. If we include normal vibrations, the 
most general form of the subordinate function will be 

F= J P(« ( X, / x,v,0,x,^). (35) 

The manner in which I have obtained equation (35) is not so direct as 
Mr. Green's method, and I have only used it for the sake of the intermediate 
equations (31), which exhibit a remarkable property of the quantities (\, p, v, 
<f>, x, ^). The direct method of deducing it is the following. Let M. Cauchy's 
elbpsoid (17) be constructed for the function V, which is homogeneous and of 
the second order with respect to (a, 0, y, u, v, w). 

Equations (16) determine the directions of molecular vibration ; in these 
equations, if (I, m, n) be substituted for (cos a, cos j3, cos y), we shall obtain 
the following: 

( Q! - E) mn + F' (n 2 - m 2 ) + Win - G'lm = 0, 
(E -F)nl + G' (I 2 - n 2 ) + F'ml - H'mn = 0, 
(F - Q') Im + IT(m 2 - l 2 ) + G'nm - F'nl = 0. 

These equations express that one of the axes of the ellipsoid is normal to the 
wave-plane, and consequently that the other two axes are contained in the wave- 
plane : by stating analytically that these conditions are true, independent of the 
position of the wave-plane, we obtain the following relations among the coeffi- 
cients of V, 

(pit) = 0, (av) = 0, (aw) = 0, 

( 7 m)=0, (7t>) = 0, (/3h>) = 0, 
(au) + 2(vw) = 0, (/3i>) + 2(uw) = 0, (yw) + 2(uv) = 0, 
(« a ) = m = (T 2 ) = 2(« 2 ) + (ft) = Hv 2 ) + (ay) = 2(v?) + (a/3). 

These fourteen equations will reduce the function V to the form 



and the Laws of Plane Waves propagated through tnem. 131 

2V= Aa> 2 + P\ + Qfi + Rv + 2F<j> + 2G% + 2H+. (35, a) 

This is the function used by Mr. Gbeen, and is a particular case of equation 
(35). The first term of this function will determine the normal vibration, and 
the last six will represent, as we have seen, transverse vibrations propagated 
according to the laws of Fresnel's wave-surface. 

If the body be homogeneous and uncrystalline, we shall have the relations, 

P=Q = R, 
F=0, 6 = 0, H=0, 

which will reduce the function V to the form 

2V = Aw i +P(\ + t i+ v ). (36) 

This function will represent homogeneous solids, liquids, and gases, and will 
be the complete function for these bodies, provided there be neither external 
forces nor pressures at the limits. If there be such forces, however, it 
will be necessary to add to the function (36), which is of the second order, 
other terms of the first order, as in equations (15). If the function (36) be 
assumed to represent all the forces engaged, the equations derived from it 
will represent the motion of a body abandoned to its own molecular actions, 
and freed from all external influence, such as gravitation, pressure of an atmos- 
phere, &c. The known properties of solids, liquids, and gases, enable us to 
determine the form of the function (36), and thus lead to the terms to be added 
in the general case for each species of body. 

It is generally admitted that a solid body, if abandoned to itself, would be 
capable of vibratory motion, and that its molecules, if displaced, would tend to 
return to their former position. A gas, if abandoned to itself, would be dissi- 
pated by the repulsive force of its molecules, so that in this case the function 
(36) should lead to an impossible result, as vibratory motion is impossible 
without the addition of pressures at the limits, or some equivalent forces. A 
liquid occupies a position intermediate between a solid and a gas ; and if we 
assume that a liquid abandoned to itself will be in a state of unstable equili- 
brium (i. e. its molecules, if displaced, will not return to their original position, 
while, if undisturbed, they will not be dissipated), we shall obtain from (36) 

s2 
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the equations of liquid motion, including friction, which have been deduced 
by various writers from different considerations. 

A liquid need not be supposed to be exactly in this state at all times ; a 
slight cohesive or a slight repulsive force may be supposed to exist among 
its molecules, according to the quantity of caloric contained in it, or other phy- 
sical circumstances, which may modify the intensity of the molecular actions. 
If such cohesive or repulsive forces be considered very small, as compared with 
the cohesive forces in a perfect solid, or the repulsive forces in a perfect gas, 
the equations deduced from the hypothesis, that these forces are zero, may still 
be used. 

We obtain from equation (36) the following: 

^=Aw-2P(a,-a), ^- = Pu, 

da ' du 

£- a. -*p („-«. £-j*. 

It is necessary and sufficient for stable equilibrium that these six forces 
should have signs contrary to the signs of (8a, 8/3, 8y, 8m, 8w, iw). If these be 
made positive, then the forces must have negative signs, and vice versa. 

Hence, for stable equilibrium, it is necessary that A and P be both negative, 
which will reduce the function (36) to the following, 

-2V = Au> 2 + P(\ + h + p). (36, a) 

Also the first three equations (changing the signs of A and P), added together, 
must be negative ; hence the condition, 

(- ZA + 4P) a, < 0. (36, b) 

4P 
If the equilibrium be stable, A cannot be less than — ; and if it be exactly 

o 

equal to this value, we shall obtain the equations pecuhar to liquids, because a 

4P 

displacement will produce no molecular force ; and if A < -7.-, a molecular force 
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will be developed which will tend to increase the displacement. The func- 
tion (36, a) will, therefore, represent solids, liquids, or gases, according as 

. AP 

A>, = ,< T . 

I shall consider, first, the equations of homogeneous solids. The function 
(36, a), substituted in equations (30), leads to the following equations of 
motion. 

e 2?- {A ~ r) fa + r W + ay + d?)' 

These equations of motion of solid bodies were first given by M. Cauchy ;* 
equations identical in form, but with the relation A = 3P between the coeffi- 
cients, had been previously obtained by M. Navier.I M. Poisson deduced 
equations identical with those of M. Navter. Mr. Green has used the equa- 
tions (37), with two independent constants, in his theory of light ;J and Mr. 
Stokes has recently called attention to the importance of retaining the two coeffi- 
cients (leaving their ratio to be determined by experiment for each solid), in a 
memoir published in the Cambridge Philosophical Society's Transactions.§ 

The relation A — ZP is a consequence of the use of definite integrals for 
the coefficients of the function V, and only represents a particular elastic solid ; 
its introduction does not alter the form of equations (37), nor does it render 
them more simple than they are in their present state. 

The additions necessary to be made to the equations of hydrodynamics, in 
order to take into account the friction of the fluid particles, have been given by 
many writers. M. Navier first stated the equations in their corrected form 
for incompressible fluids.|| M. Poisson has treated of the subject in a me- 

* Exercices des Mathematiques, torn. iii. p. 180. 

f Memoires de l'Institut, torn. vii. p. 389- 

\ Transactions of Cambridge Philosophical Society, torn. vii. p. 1 1 . 

§ Vol. viii. part 3. || Memoires de l'Institut, torn. vi. p. 414. 
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moir published in the Journal de VEcole Polytecknique ;* and, more recently, 
M. Babbe de Saint Tenant) - and Mr. Stokes J have "written on the friction of 
fluids in motion. 

The quantities to be added to the ordinary equations of hydrodynamics are 

the right hand members of equations (37), introducing the relation I A = -^- J , 

which expresses that a displacement produces no molecular force. 

The equations of the motion of gases have been already given, on the sup- 
position that there is no tangential action, which is equivalent to assuming that 
the vibrations are normal, and therefore (P = 0). If we wish to take account 
of the friction of gases, we should use the equations which have just been indi- 
cated for liquid motion. 

section v. compaeison of mechanical theoeies of light. 

It is well known that different mechanical theories have been proposed to 
account for the phenomena of the movement of light in crystalline bodies, and 
that, although these theories differ in their fundamental hypotheses, yet, to some 
extent, they agree in representing the most obvious phenomena of double re- 
fraction. The laws of wave-propagation in crystals, are geometrical conse- 
quences of the properties of Feesnel's wave-surface; and no mechanical theory 
of light can be considered as even an approximation to the truth, unless it con- 
tains, as a deduction from its hypotheses, the wave-surface of Fbesnel. But 
it would be an error to conclude that any theory is correct, which satisfies this 
condition. There are, in fact, three different theories which satisfy this funda- 
mental condition, and it is evident that they cannot all be true. The first of 
these theories was propounded by Fbesnel himself, in his memoir on Double 
Refraction.§ It is based on the hypothesis, that the luminiferous ether is com- 
posed of attracting and repelling molecules. The form of wave-surface known 
as Feesnel's is deduced by its author from this hypothesis, with the peculiarity 
that the vibrations of the molecules are perpendicular to the plane of polarization. 

* Cahier xx. p. 139. t Comptes Eendus, torn. xvii. p. 1240. 

J Transactions of Cambridge Philosophical Society, vol. viii. part 3. 
§ Memoires de l'Institut, torn, vii, p. 45. 
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M. Cauchy afterwards gave the general equations peculiar to such a system, 
and deduced from them Fbesnel's wave-surface, as a first approximation to what 
he considered as the more accurate laws of wave-propagation.* 

In the year 1839, Mr. Gbeen presented to the Cambridge Philosophical 
Society a memoir, in which, by a modification of M. Cauchy's equations, he 
obtained Fbesnel's wave-surface as an exact deduction from the theory, f 
This modification consists, as I have already stated, in restricting the system to 
propagate normal and transverse vibrations. In M. Cauchy's or Mr. Green's 
theory, the vibration of the molecules is parallel to the plane of polarization. 
In the same year Professor Mac Cullagh presented to this Academy a mecha- 
nical theory of light, not founded on the hypothesis of attracting and repelling 
molecules. The vibrations in this theory also are parallel to the plane of pola- 
rization, and the form of the wave-surface is that given by Feesnel. These 
three theories of fight, therefore, agree, so far as the laws of wave-propagation 
are concerned ; and, excluding Fresnel's theory from the comparison (as the 
vibrations perpendicular to the plane of polarization make it distinct from the 
other two theories), there remain the mechanical theories of Mr. Gbeen and Pro- 
fessor Mac Cullagh, which are identical so far as the laws of wave propagation 
are concerned. The two theories are, however, really different in their funda- 
mental assumptions ; and this remarkable agreement in the laws of wave-propa- 
gation deduced from them admits of a simple explanation. I propose to account 
for the agreement, and to suggest the direction in which we should look for a 
true experimentum cruets between them. 

The function V used by Mr. Gbeen, when reduced to its simplest form, 
will be 

-2jr=Aw* + P\+Qn + Bv, (38) 

and the simplest form of Professor Mac Cullagh's equations will be derived 
from the function 

-2V=PX*+QY 2 + RZ\ (39) 

It is evident from what I have stated in the first section, that (\, ft, v) will 

* Memoires de l'Institut. torn, x., 1830. 

f Transactions of Cambridge Philosophical Society, vol. vii. p. 121. 



136 The Rev. S. Hatjghton on a Classification of Elastic Media, 

produce the same terms in the differential equations of motion, as X 2 , Y 2 , Z 2 ; 
for the squares will be the same in each, and the rectangles will be 

(2/3372-4^73, 27,03-40,73, 20JJ8, - 4o,j8 2 ), and (-2&7 S , -27,03, -2^). 

These two sets of rectangles will produce the same terms in the equations of mo- 
tion, since we may transpose the differentiations without affecting the result, so 
far as the laws of propagation are concerned. The surfaces of wave-slowness de- 
duced from (38) and (39) may be obtained immediately from equation (18). 
Equation (38) will give the following: 

P = Al 2 + #m 2 +Qn 2 ; F'=(A-P)mn; 

Q! = Am 2 + Pn 2 + Rl 2 ; G' = (A -Q)ln; (40) 

R = An 2 + Ql 2 + Pm 2 ; H T ={A-R) Ivn. 

And similarly from equations (39) will be found 

P = Rm 2 + Qn 2 ; F' = - Pmn ; 

Q! = Pn 2 + Rl 2 ; G' = - Qln ; (41) 

R=Ql 2 +Pm 2 ; R^-Rlm. 

These equations differ from the former only by not containing .4. 
The equation of wave-slowness (18) derived from (40) is 



\A(x 2 + y I +z 2 )-l\ 
x { (^y 2 +z 2 )(QRic 2 +PRy 2 +PQz 2 )-(Q+R)^-(P+R)y 2 -(P+Q)z 2 +l }=0. 



(42) 



The first factor of this equation represents a sphere whose radius is —t-t, and 

belongs to the normal vibration; the second factor is the equation of Fbesnel's 
wave-surface, and in it the vibrations are transversal. The equation of wave- 
slowness deduced from Professor Mac Cuixagh's function (39) will be the last 
factor of (42). So far, therefore, as the laws of wave-propagation are concerned, 
the functions (38) and (39) are equivalent, with this difference, that the func- 
tion (38) introduces a normal wave, which does not enter into the equations 
derived from (39). It might be thought at first sight that we are at liberty to 
make A = 0, and thus reduce the function (38) to a function representing 
nothing but transverse vibrations ; this, however, cannot be admitted, for as 
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(38) represents a body whose molecules act in the line joining them, a wave of 
normal compression is always possible. This will be rendered more evident by 
considering the conditions at the limits. 

Let the limiting surface separating two bodies be the plane (a?, y), then the 
equations of condition (11) will become, for the function (38), 

GV. = QV.'; £ = £'; 

Pu' = P'u' o '; % = %; ( 4S ) 

A'w' a + Q< + Pfr = A"«o'+ Q'tf + P'fi'i to = £'• 

These equations are equal in number to the unknown quantities, provided nor- 
mal waves be included, because the unknown quantities of the problem are the 
intensities of the reflected and refracted waves ; it is impossible, therefore, for 
exclusively transverse waves to be produced by reflexion or refraction in such 
a body as (38) defines : in order to obtain unknown quantities whose number 
shall be equal that of the necessary conditions of the mechanical problem, we must 
introduce normal vibrations. The conditions at the limits deduced from ( 3 9 ) are 

*°l -» = V ■* > So = hO » 

P'X'^P'X','- V ' = V ' '. (44) 

The additional hypothesis made by Professor Mac Cullagh, that the density 
of the medium is the same in the two bodies, reduces these equations to four. 
Accordingly, with this hypothesis, there is no necessity to have recourse to normal 
waves, as there will be four intensities to be determined in the transverse waves. 

From these considerations it appears, that the experimenta cruris between 
the rival theories of light must be sought for among the laws of reflexion and 
refraction ; but unfortunately these laws are not known with sufficient accuracy 
to enable us to decide the question. Mr. Green's theory contains the common 
laws of reflexion at the surfaces of ordinary media as first approximations, while 
Professor Mac Cuixagh's system has the advantage of giving these laws as exact 
results ; nothing, however, but more accurate experiments can decide whether 
the approximation or the exact result be most in accordance with the truth ; and 
as these experiments involve considerations of the intensity of light, it would be 

VOL. XXII. T 
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difficult to make them with sufficient accuracy. The present state of the wave 
theory of light certainly suggests grave doubts as to the nature of the foundation 
on which the whole system is based. We first assume the existence of an un- 
known medium, whose existence must remain unproved and unprovable by us ; 
then, from supposed properties of this unknown medium, we deduce the laws of 
propagation, &e. Here a new difficulty arises; for we find several different theories 
capable of explaining the laws of propagation, and explaining with more or less 
exactness the most obvious of the laws of reflexion and refraction. How are we 
to decide among these conflicting theories ? Are we to assume, with M. Caucht, 
that the observed laws of polarized light occupy, with respect to the mathematical 
laws deduced from his theory, the same position that the laws of Kefleb stand in 
with respect to the more accurate laws of planetary motion? or are we to assume 
that theory to be correct which agrees accurately with the common formulas 
for reflected light, when it is well known that these formulas themselves are 
doubtful for highly refracting substances? It appears certain, that we do 
not yet possess experimental knowledge sufficient to enable us to determine 
which of the theories of light is correct, or whether any of them be so. In a 
general point of view Professor Mac Cullagh's theory possesses an important 
advantage, as compared with other theories. It contains no inexplicable normal 
wave, and does not render this difficult subject still more intricate, by the in- 
troduction of a useless vibration. It is greatly to be desired, that the attention 
of experimentalists were directed to the necessity which exists for more accurate 
and general researches into the laws of crystalline reflexion and refraction, and 
that the surface of Fresnel were placed upon a purely experimental basis. From 
such researches, carefully conducted, might be deduced the geometrical laws of 
double refraction, and a foundation be laid for a complete and positive theory 
of the laws of polarized light. 



